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Abstract 

This manuscript is devoted to investigate Bianchi Type I universe 
in the context of f{R,T) gravity. For this purpose, we explore the 
exact solutions of locally rotationally symmetric Bianchi type I space- 
time. The modified field equations are solved by assuming expansion 
scalar 0 proportional to shear scalar a which gives A = B^, where 
A, B are the metric coefficients and n is an arbitrary constant. In 
particular, three solutions have been found and physical quantities are 
calculated in each case. 
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1 Introduction 

Recent observations from astrophysical data have unfolded an amazing pic¬ 
ture of expanding universe. The cosmic acceleration is well supported by high 
red-shift supernovae, cosmic microwave background anisotropy and galaxy 
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clustering [T]. Universe seems to be filled with exotic cosmic fiuid known as 
dark energy having strong negative pressure. It constitutes almost 70 % of 
the total energy budget of our universe. We can describe dark energy with 
an equation of state (EoS) parameter u = p/p, where p and p represent the 
energy density and pressure of dark energy. It has been proved that the ex¬ 
pansion of the universe is accelerating when tc — 1 [2]. The phantom like 
dark energy is found to be in the region where a; < — 1. The universe with 
phantom dark energy ends up with a finite time future singularity known as 
cosmic doomsday or big rip j3]. The modified theories of gravity seem attrac¬ 
tive to explain the phenomenon of dark energy and late time acceleration. It 
is now expected that the issues of cosmic acceleration and quintessence can 
be addressed using higher order theories of gravity |1]. 

Many generalizations of Einstein field equations have been proposed in 
the last few decades. /(T) theory of gravity is an example which has been 
recently developed. This theory is a generalized version of teleparallel gravity 
in which Weitzenbock connection is used instead of Levi-Civita connection. 
The theory seems interesting as it may explain the current cosmic accel¬ 
eration without involving the dark energy. Some researchers have done a 
considerable amount of work in this theory so far [5]. Another extended 
theory known as f{R) theory of gravity has also attracted attention of the 
researchers in recent years. The f{R) theory is actually an extension of stan¬ 
dard Einstein Hilbert action involving a function of the Ricci scalar R. The 
cosmic acceleration may be justified by involving the term 1/R which is re¬ 
quired at small curvatures. The f{R) theory seems to be most appropriate 
due to important f{R) models in cosmological contexts. 

Some viable f{R) gravity models [6] have been suggested which justify 
the unification of early-time inflation and late-time acceleration. The dark 
matter problems can also be addressed using viable f{R) gravity models 
Starobinsky [8] gave first complete inflationary R + R^ model viable with the 
observational data. The cosmological model without additional singularities 
can be constructed using f{R) gravity in which both inflation in the early 
universe and dark energy in the present universe is described [H]. Hendi 
and Momeni m investigated black hole solutions in f{R) theory of gravity 
with conformal anomaly. Jamil et ah to explored f{R) tachyon model using 
Noether symmetries Capozziello et ah |T2] have proved that dust matter 
and dark energy phases can be achieved by finding the exact solutions using 
a power law f{R) cosmological model. A reasonable amount of work has 
been done so far in this theory |T3]. Some review articles [2] may be useful 
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to have a better understanding of the theory. 

Recently, Harko et ah [15] proposed a new generalizations known as 
f{R, T) theory of gravity in which R is the scalar curvature and T denotes the 
trace of the energy-momentum tensor. Jamil et ah [16] reconstructed some 
cosmological models in f{R, T) gravity where it was proved that the dust fluid 
reproduced ACDM. Sharif and Zubair im gave the reconstruction and stabil¬ 
ity conditions of f{R,T) gravity with Ricci and modihed Ricci dark energy. 
The same authors [18] discussed the laws of thermodynamics in this theory. 
However, it has been established that the hrst law of black bole thermody¬ 
namics is violated for f{R,T) gravity [19]. Santos [20] investigated Godel 
type universe in the context of f{R, T) modihed theories of gravity. Houndjo 
fIT\ reconstructed f{R,T) gravity by taking f{R,T) = fi{R) + f 2 {T) and it 
was shown that f{R, T) gravity allowed transition of matter from dominated 
phase to an acceleration phase. Harko and Lake [22] investigated cylindrically 
symmetric interior string like solutions in f{R,Lm) theory of gravity. In a 
recent paper [23], we explored the exact solutions of cylindrically symmetric 
spacetime in f{R, T) gravity and recovered two solutions which corresponded 
to an exterior metric of cosmic string and a non-null electromagnetic held. 

Isotropization is an important issue to discuss whether the universe can 
result in isotropic solutions without the need of hne tuning the model param¬ 
eters [23]. The universe seems to have an isotropic and homogeneous geom¬ 
etry at the end of the inhationary era [23] • However, the class of anisotropic 
geometries has gained popularity under the light of the recently announced 
Planck Probe results [26] . Further, it is believed that the early universe may 
not have been exactly uniform. Therefore, inhomogeneous and anisotropic 
models of universe have an important role in theoretical cosmology. This 
prediction motivates us to describe the early stages of the universe with the 
models having anisotropic background. Thus, the existence of anisotropy 
in early phases of the universe is an interesting phenomenon to investigate. 
Bianchi type models are among the simplest models with anisotropic back¬ 
ground. Many authors [2Z] investigated Bianchi type spacetimes in different 
contexts. Kumar and Singh [2H] solved the field equations in the presence 
of perfect fluid using Bianchi type I spacetime in general relativity (GR). 
Moussiaux et ah [29] investigated the exact solution for vacuum Bianchi 
type III model in the presence of cosmological constant. Bianchi type III 
string cosmology with bulk viscosity has been studied by Xing-Xiang [30] . 
He assumed expansion scalar proportional to the shear scalar to hnd the so¬ 
lutions. Wang [21] explored string cosmological models in Kantowski-Sachs 
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spacetime. Magnetized Bianchi type III massive string cosmological models 
in GR have been investigated by Upadhaya [52]. Hellaby [33] gave a review 
of some recent developments in inhomogeneous models and it was concluded 
that the universe is inhomogeneous on many scales. 

The investigation of Bianchi type models in modihed or alternative the¬ 
ories of gravity is also another interesting topic of discussion. Perfect fluid 
solutions using Bianchi type I spacetime in scalar tensor theory have been 
explored by Kumar and Singh [33] • Singh et ah [33] studied Bianchi type 
III cosmological models in scalar tensor theory. Adhav et ah [36] found 
an exact solution of the vacuum Brans-Dicke held equations for a spatially 
homogeneous and anisotropic model. FRW cosmologies in f{R) gravity have 
been investigated by Paul et ah [3Z]. Bianchi type / model in f{R) grav¬ 
ity was studied where it was shown how to integrate anisotropic degrees of 
freedom explicitly and to reduce the problem to one differential equation 
for the volume factor [33]. Vacuum and non-vacuum solutions of Bianchi 
types I and V spacetimes in metric f{R) gravity have been explored [391130] . 
Sharif and Kausar [IT] investigated non-vacuum solutions of Bianchi type 
VI universe by considering the isotropic and anisotropic huids as the source 
of dark matter and energy. In a recent paper, we have explored Bianchi 
type / cosmology in f{R,T) gravity with some interesting results [12]. It 
was concluded that equation of state parameter w—)■—last—)-oo which 
suggested an accelerated expansion of the universe. Thus it is hoped that 
f{R,T) gravity may explain the resent phase of cosmic acceleration of our 
universe. This theory can be used to explore many issues and may provide 
some satisfactory results. 

In this paper, we explore the exact solutions of locally rotationally sym¬ 
metric (LRS) Bianchi type / spacetime in f{R,T) gravity. The held equa¬ 
tions are solved by assuming expansion scalar 6 proportional to shear scalar 
a which gives A = where A, B are the metric coefficients and n is an ar¬ 
bitrary constant. The plan is planned as follows: Field equations in f{R,T) 
gravity are briehy introduced in section 2. In section 3, the solutions of the 
held equations are investigated along with some important physical param¬ 
eters. Last section is used to summarize and conclude the results. 
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2 f{R, T) Gravity Formalism 


The action for f{R,T) theory of gravity is given by 


S= y/^{—f{R,T) + Ljd\ 


( 1 ) 


where g is the determinant of the metric tensor and Lm is the usual 
matter Lagrangian. It would be worthwhile to mention that if we replace 
f{R,T) with f{R), we get the action for f{R) gravity and replacement of 
f{R, T) with R leads to the action of GR. The f{R, T) gravity field equations 
are obtained by varying the action S in Eq.([T]) with respect to the metric 
tensor 


fR{R,T)R,,--f{R,T)g^,-{V^V,-g^,n)MR,T) = r)(r^,+0^,), 


where denotes the covariant derivative and 


( 2 ) 


P) _ a/3 


dR ’ ^ ffr 

Contraction of Eq.(l2]) yields 

/«(R, T)R + 3n/^(R, T) - 2f{R, T) = kT - MR, T){T + 0), (3) 


where 0 = 0^^. This is an important equation because it provides a rela¬ 
tionship between Ricci scalar R and the trace T of energy momentum tensor. 
Using matter Lagrangian the standard matter energy-momentum tensor 
is derived as 

M = (p + p)Uf,u^ - pg^y, (4) 


satisfying the EoS 


p = wp, 


(5) 


where 0, 0, 0) is the four-velocity in co-moving coordinates and 

p and p denote energy density and pressure of the fluid respectively. Perfect 
fluids problems involving energy density and pressure are not any easy task 
to deal with. Moreover, there does not exist any unique dehnition for matter 
Lagrangian. Thus we can assume the matter Lagrangian as = —p which 
gives 

Qfiu = -pPfiu - Mu, ( 6 ) 
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and consequently the field equations ([2]) take the form 


(7) 

It is mentioned here that these field equations depend on the physical nature 
of matter field. Many theoretical models corresponding to different matter 
contributions for f{R,T) gravity are possible. However, Harko et ah [15] 
gave three classes of these models 

r R+2fiT), 

f{R,T)={ f,{R) + MT), 

[ U{R) + UR)UT). 

In this paper, we consider the first and second class only to explore the exact 
LRS Bianchi I solutions. 


3 Exact LRS Bianchi Type / Solutions 

Here we first develop some important cosmological parameters and field equa¬ 
tions for LRS Bianchi type I spacetime and then find the exact solutions of 
field equations for constant and non-constant curvature case. 


3.1 LRS Bianchi Type I Spacetime 

The line element of LRS Bianchi type I spacetime is given by 

= dt^ — A^{t)dx^ — B'^{t)[dy‘^ + dz^]i (8) 


where A and B are cosmic scale factors. The corresponding Ricci scalar turns 
out to be 


R 


,A 2B 2AB 

+ "b + ^ + fij'- 


(9) 


where dot denotes derivative with respect to t. The average scale factor a 
and the volume scale factor V are defined as 


a = ^AB^, V = a^ = AB^. 


( 10 ) 
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The average Hubble parameter H is given in the form 


H 




The expansion scalar 9 and shear scalar a are defined as follows 


(11) 


where 


9 = 


a = 


„ A B 

= I + 


1 




-[- - -] 




( 12 ) 

(13) 

(14) 


hfiu = gfiu — u^Uiy is the projection tensor. Now we explore the solutions of 
the held equations for two classes of f{R,T) models. 


3.2 f{R,T) = R + 2f{T) 

For the model f{R, T) = R + 2/(T), the held equations become 


R 


fiiy 



Rl'fiu + 2 /t(T)T^j, + 


/(T) + 2p/t(T) 


Qrv 


(15) 


Here we hnd the most basic possible solution of this theory due to the com¬ 
plicated nature of held equations. For the sake of simplicity, we use natural 
system of units (G = c = 1) and /(T) = AT, where A is an arbitrary con¬ 
stant. In this case, the gravitational held equations take the form similar to 
GR 

R/ii. - ]^R9i,u - A(T -h 2p)g^^ = (Svr -h 2X)T^^, (16) 

where the term A(T -|- 2p) may play the role of cosmological constant A of 
the GR held equations. It would be worthwhile to mention here that the 
dependence of cosmological constant A on the trace of energy momentum 
tensor T has already been proposed by Poplawski [43] and the cosmological 
constant in the gravitational Lagrangian is a function of T. Gonsequently 
the model was named as “A(T) gravity”. It has been proved that recent 
astrophysical data favor a variable cosmological constant which is consistent 
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with A(T) gravity. A(T) gravity has been shown to be more general than the 
Palatini f{R) gravity |11]. Now using Eq. fllbp . we obtain a set of differential 
equations for LRS Bianchi type / spacetime, 

(Svr + 3A)p - Ap, (17) 

(Svr + 3A)p — \p, (18) 

(Svr + 3A)p — \p. (19) 

Thus we have three differential equations with four unknowns namely A, B, 
p and p. Adding Eqs. flTT)) and flT 8 |) gives 


2AB B^ 
~AB^^ 
2B 

A B AB 
A ~ B ~ AB 


2AB 2B 
AB ~ 


(87r + 2A)(p + p). 


Similarly, addition of Eqs. flTTl) and flT^ yields 


AB A B B"^ 
AB “ A “ B ^ 


(87r + 2A)(p + p). 


Subtracting Eqs. fl2UD and ([21]), it follows that 


( 20 ) 


( 21 ) 


AB A B B'^ 

'AB^ ~A~ ~B~ 


( 22 ) 


Now we are left with only one differential equation and two unknowns. There¬ 
fore, we need an additional constraints to Eq. fl22|) . Here we use a physical 
condition that expansion scalar 9 is proportional to shear scalar a which 
provides 

A = B", (23) 

where n is an arbitrary real number and we consider n 7 ^ 0,1 for non-trivial 
solutions. The physical reason for this assumption is justihed as the observa¬ 
tions of the velocity red-shift relation for extragalactic sources suggest that 
Hubble expansion of the universe may achieve isotropy when | is constant 
[U]. Collins |1B] provided the physical signihcance of this condition for per¬ 
fect fluid with barotropic EoS. In literature ra, Bii-iMi, many authors have 
proposed this condition to hnd the exact solutions of held equations. 







(24) 


Thus using Eq. fE^ . Eq. (12^ take the form 


p> tj2 

^ + (" + 1 )^ = 0 . 


which yields a solution 


5 = ci[(n + 2)t + C2]~+2, (25) 

where ci and C 2 are constants of integration. Thus, the solution metric takes 
the form 

ds^ = dt^ — ci^'^[{n + 2)t + C2)\'^dx^ — Ci^[{n + 2)t + C2)\^[dy'^ + dz^]. (26) 
The volume scale factor turn out to be 


V = a^ = c^^+\n + 2)t + C2]. 


(27) 


The expansion scalar and the shear scalar become 


e 


n + 2 

{n + 2)t + C 2 


1 n — 1 
3 _ (77. + 2)t + C 2 


(28) 
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It is mentioned here that the isotropy condition, i.e., ^—)■ 0 as t —)■ cx3, is 
satished in this case. It can be observed from Eqs. fl27j) and fl28l) that the 
spatial volume is zero at t = 0 while the expansion scalar is inhnite, which 
suggests that the universe starts evolving with zero volume at t = 0, i.e. big 
bang scenario. It is further observed that the average scale factor is zero at 
the initial epoch t = 0 and hence the model has a point type singularity [52] . 
The energy density and pressure of the universe take the form 


(1 + 2n) 

2(A + 47r)[(?7, + 2)f + 02 ]^ 


(29) 


which suggest that equation of state parameter a; = 1 corresponding to stiff 
fluid universe. The average Hubble parameter turn out to be 


Therefore 


H 

3[(n. + 2)t + C 2 ] 

(30) 

H (71 + 2 )fg + C2 

Hq {ji + 2 )f + C2 

(31) 
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where Hq is the present value of Hubble’s parameter. The redshift for a 
distant source is directly related to the scale factor of the universe at the 
time when the photons were emitted from the source. The scale factor a and 
redshift 2 ; are related through the equation 


a = 


ao 

1 + z’ 


(32) 


where oq is the present value of the scale factor. Thus we obtain 

(n + 2)to + C 2 ® 

(n + 2)f + C 2 


Oo 


= 1 + 2 = 


(33) 


Using Eqs. (1^ and fl33H . we obtain the value of Hubbles parameter in terms 
of redshift parameter 

H = Ho(l + zf. (34) 

According to the Hubble Law, the distance of a given galaxy is proportional 
to the recessional velocity as measured by the doppler red shift. Thus, the 
value of Hubbles parameter in terms of redshift has much importance in 
astrophysical contexts. 

The deceleration parameter q in cosmology is the measure of the cosmic 
acceleration of the universe expansion and is dehned as 


da 

It is mentioned here that the behavior of the universe models depend upon 
the sign of q. The positive deceleration parameter provides a decelerating 
model while the negative value corresponds to inflation. For this solution, 
the value of deceleration parameter turns out to be g = 2 which suggests a 
decelerating model of universe. 

Universe models closed to ACDM can be described using the cosmic jerk 
parameter j, a dimensionless third derivative of the scale factor with respect 
to the cosmic time [53]. The value of jerk parameter is constant for flat 
ACDM model. The jerk parameter is defined as 



1 d 
^ ” IPa' 


(36) 


The expression for jerk parameter in terms of deceleration parameter turns 
out to be 


j =q + 2q^ - 


(37) 
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Thus we obtain j = 10 in the case of our solution. It would be worthwhile 
to mention here that this solution gives R = 0 for n = — |. In this case, the 
solution metric takes the form 

ds'^ = dt^ - [ci(^t + C 2 )\~^dx^ - [ci(^t + C 2 )]^[dy^ + dz^]. (38) 


Without loss of generality, we take C 2 = 0 and re-define the parameters, i.e., 

^ y and 


2x^ 
3ci3/2 

form 


3/ 9ci3^3 

4 


_3 / 9ci323 


ds^ = dt^ — t 3dx'^ — t3[dy‘^ + dz'^), 
which is exactly the same as the well-known Kasner’s metric 


z, the above metric takes the 

(39) 


Now we discuss the possibility of solutions for a non-linear forms of f{T). 
We assume /(T) = AT^ so that Eq. flT^ takes the form 

Rfiu - = (Stt 4AT)T^j. AT(T 4:p)gf,^,. (40) 

Now using Eq. fl4U]) . we obtain a set of independent differential equations for 

= Svrp -I- 5Ap^ — 14App — 3Ap^, (41) 

= Sirp — + 6Xpp — Xp^, (42) 

= Sup — 9Xp^ + 6Xpp — Xp^. (43) 

These equations yield the same solution metric as given by Eq. fl38ll . How¬ 
ever, in this case the energy density can be obtained by solving the equation 

2 ^ ^ 

^ ^ 2A(1 -3a;)^~ 

The quadratic equation is due 




8A(1 — 3a;) [(u -|- 2)t -|- 02 ]^ 
to non-linear form of /(T). 


= 0 . 


(44) 


LRS Bianchi type / spacetime, 

2AB 
~AB ^ 

2B B^ 

~ B^ 

A B AB 
H “ AB 
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3.3 f{R,T) = MR) + f2{T) 

Now we explore the solutions with more general class. Here the held equa¬ 
tions for the model f{R,T) = fi{R) + f 2 {T) become 


= [«+/2T(T)]T^.+[p/2r(T)+i/2(T)]^7^.. 

(45) 

Contracting the held equations fl45|l . we obtain 

RhniR) - 2MR) + 3Dh^{R) = nT + 2/2(T) + [T + 4p]/2^(T). (46) 

Using this, we can write 

^ 3D/i^(i?) + RhniR) -^T- 2/2(T) - [T + Ap]hAT) _ 


Inserting this in Eq. (l45ll . we get 

fin{R)RRu - V,VJin{R) -{n + / 2 t(T))T^. 
9r,u 

RfiRiR) - □/!«(«) - Kr - Tf2r(T) 


(48) 


Since the metric (jS]) depends only on t, one can view Eq. fHHj) as the set 
of diherential equations for fijiit), / 2 'r(^); B) P and p. It follows from 
Eq. fHSl) that the combination 


, _ - (k + hT(T))T„, 

Qrr, 


(49) 


is independent of the index p and hence — Ay = {] for all p and v. Thus 
Ho — Hi = 0 yields 


2B 2AB AfUR) /UR) 
B AB ^ AR„(R) /,„(fl) 


' k +um 

. Rr(R) 


(p + p) = 0. 


(50) 


Also, Hq — H 2 = 0 provides 


A B AB B^ BR„(B.) fUn) 
B^ AB^ B^- BR^(R) Rg{R) 


^ + / 2 t(^) 

. finiR) . 


{p + p) (51) 
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Now we have two differential equations with six unknowns namely A, B, 
/i(-R), / 2 (r), p and p. Here we also use A = so that Eqs. (ISUD and (IHTD 
take the form 


2B 


B^ + Bf,^{R) 


B 


B^ 


in 




BfiniR) , hniR) 


B 


+ 


Bf,j,iR) 


+ 


^ + / 2 t(^) 

. MR) . 

^ + /2 t(^) 

. fiAR) . 


(p + p)= 0,(52) 
(P + P) = 0.(53) 


Case I: Exponential Law Solutions 


It has been proved that dark matter and dark energy phases can be achieved 
by Ending the exact solutions using a power law f{R) model [ 12 ]. So it 
would be interesting to assume fi{R) in power law form to solve the field 
equations. We follow the approach of Nojiri and Odintsov [55] and take the 
assumption /i^(i?) oc foR^, where /o is an arbitrary constant. So in this 
case the addition of Eqs. (15211531) yields 

(n^ + 2n^ + 2n + 2)B^ + (2n^ + 4n + 3)BB^B + (n + 2)B^B^ - 
2m{n^ + n + 1).B^ + m{2in? + n)BB^B + m{n + 2)B^BB = 0. (54) 

It is interesting to notice that this equation admits an exponential solution 
of the form 

H(^) = e"^*+"^ (55) 

where C 3 and C 4 are arbitrary constants. The exponential solution is satisfied 
with the constraint equation 


+ ArR -|- 7n -|- 6 = 0. 

( 66 ) 

The solutions of this equation turn out to be 

n = — 2 , —1 ± i\/ 2 . 

(57) 


It is mentioned here that the real value of n gives a constant Ricci scalar 
while we obtain a non-constant Ricci scalar for the complex values of n. We 
discard the imaginary case and consider the real value of n to get a physical 
solution 

ds^ = dR - + dz^). (58) 
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The average Hubble parameter turn out to be zero here. All other dynamical 
quantities like volume scale factor of universe, expansion scalar 6 and shear 
scalar a are constant for this solution. Energy density and pressure of the 
universe are related by the equation 


P + P 


(_6)™+7oC32 

« + /2Tm 


(59) 


Many expressions for pressure and energy density can be evaluated for dif¬ 
ferent choices of / 2 (T). For example when / 2 (T) = AT^ and using Eq.([5]), we 
obtain 


2 ^ (-6)"^+VoC3^ _ 

^ 2A(l-3n;)^ 2A(1 + n;)(l - Sw) 


(60) 


which is quadratic in p and one can work out its roots to get energy density. 


Case II: Power Law Solutions 


Here we assume that the solution is in power law form, i.e. B{t) = (cst-t-Cg)^, 
where C 5 , Cg and k are arbitrary real constants with k ^ 0. Using Eq. dSTj) . 
we obtain a constraint equation 

{n^ + An^ + 7n + 6)k'^ — [2m{n‘^ + 2n + 3)+ 271^ + 6n + 7]k + {l + 2m){n + 2) = 0. 

(61) 

This equation is important because it will be used to reconstruct different 
forms of fi{R) models with suitable solutions of field equations. For example, 
here we investigate the solution for n = —1. In this case, Eq. fl 6 ip gives 

k = 2m -|- 1 , m 7 ^ 0 , 1 . 

So the solution metric takes the form 

ds^ = dt^ — {c^t + dx‘^ — {c^t + + dz^). 

The volume scale factor and average Hubble parameter become here 

r = 0= = (c,t + H = 

3(C5t Cg) 

The expansion scalar and shear scalar turn out to be 

^ _ C5(2m -F 1) ^2 _ 4. C5(2m -F 1) ^ 

Cgt -|- Cg 3 Cgt -|- Cg 


(62) 

(63) 

(64) 

(65) 
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The isotropy condition ^—)■ 0 as t —)■ oo, is also satisfied in this case. Using 
Eg.lpTl). we get 

H _ c^to + Ce Oo 

Hq C^t + Cg Oi 


Cgto + Cg 
Cgt + Cg 


( 66 ) 


Thus the value of Hubbles parameter in terms of redshift parameter turns 
out to be 


H = Hq{1 + 2;)2(2 ™+i) . 
The deceleration parameter in this case becomes 

1 — 4m 


2 (2m + 1) 


(67) 

( 68 ) 


while the jerk parameter is given by 

8(4m^ — 5m + 1) 

^ ( 2 m+ 1 )^ 


(69) 


By observing Eqs. (167116^ . it is clear that singularity occurs at m = — 
Further, the Ricci scalar turns out to be 


2 c5^(1 + 2m) (1 + 4m) 

(cgt + Cg)2 


(70) 


For m = —I or m = —the Ricci scalar turns out to be constant, i.e. i? = 0. 
m = — i corresponds to Minkwoski spacetime while for m = — | the metric 
takes the form 


ds^ = dt^ — {c^t + Cg) ^dx^ — {c^t + CQ){dy‘^ + dz^). (71) 


This solution gives a point singularity at t The Ricci scalar remains 

non-constant for m ^ For a special case when m = —1, /i(R) takes 

the logarithmic form 

/i(R)oc/oln|R|+C7, (72) 

where cy is an integration constant. Here energy density and pressure of the 
universe are related by the equation 


P + P 


fo 

^ + / 2 t(^) 


(73) 


Here we can also calculate many expressions for energy density depending 
upon the value of / 2 (T). 
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4 Concluding Remarks 

This paper is devoted to study Bianchi type cosmology in f{R,T) gravity. 
We explore the exact solutions of field equations for LRS Bianchi type I 
spacetime. Since the held equations are highly nonlinear and complicated, 
we use the assumption that the expansion scalar 6 is proportional to the 
shear scalar a to solve them. It gives A = where A, B are the metric 
coefficients and n is an arbitrary constant. Mainly we have explored three 
solutions of modihed held equations using diherent assumptions. 

The hrst solution is obtained for the model /(i?,T) = R + 2/(T). The 

2 

isotropy condition, i.e. ^—)■ 0 as t —>■ oo, is satished for the solution. The 
spatial volume is zero at t = 0 and the expansion scalar is inhnite, which 
suggests that the universe starts evolving with zero volume at t = 0, i.e. 
big bang scenario. The average scale factor turns out to be zero at the 
initial epoch t = 0 and hence the model has a point type singularity [52] . 
The expressions for energy density and pressure suggest that EoS parameter 
uj = 1 which corresponds to stih huid universe. The deceleration parameter 
q turns out to be g = 2 which suggests a decelerating model of universe. We 
have also calculated jerk parameter j = 10 in the case of our solution. It is 
worth mentioning here that this solution gives i? = 0 for n = — | corresponds 
to well-known Kasner’s solution already available in GR [M] . 

We have also explored the more general solutions of held equation by 
considering the model /(R, T) = fi{R) + f 2 (T). Moreover, we have not used 
any conventional assumption like constant deceleration parameter or varia¬ 
tion law of Hubble’s parameter to investigate the solutions in this case. In 
particular, exponential law and power law solutions have been investigated 
for this model. Foe exponential law case, the average Hubble parameter turn 
out to be zero and all other dynamical quantities like volume scale factor of 
universe, expansion scalar 6 and shear scalar a are constant. It would be 
worthwhile to mention here that when f 2 {T) = 0, this class corresponds to 
f{R) gravity model. For / 2 (T) ^ 0, different expressions for energy den¬ 
sity can be generated with different choices of f 2 {T) models. Power law 
solution provides a non-constant scalar curvature and thus many important 
f{R) models can be reconstructed. As a special case, we have developed an 
important logarithmic f{R) model. 
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